Entanglement production due to quench dynamics of an anisotropic XY chain in a 

transverse field 
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We compute concurrence and negativity as measures of two-spin entanglement generated by a 
power-law quench (characterized by a rate and an exponent a) which takes an anisotropic 
XY chain in a transverse field through a quantum critical point (QCP). We show that only spins 
separated by an even number of lattice spacings get entangled in such a process. Moreover, there 
is a critical rate of quench, r~^, above which no two-spin entanglement is generated; the entire 
entanglement is multipartite. The ratio of the entanglements between consecutive even neighbors can 
be tuned by changing the quench rate. We also show that for large r, the concurrence (negativity) 
scales as \fajT [a/r), and we relate this scaling behavior to defect production by the quench 
through a QCP. 

PACS numbers: 03.67.Mn, 73.43.Nq, 64.60.Ht, 75.10.Jm 



I. INTRODUCTION 

The role of entanglement in the theory of quantum 
phase transitions has been a subject of recent studies 0]. 
A number of such works have pointed out that entan- 
glement can be used as a tool to characterize quantum 
phase transitions for both clean and disordered systems 
ISs3,s^i5ii6] . These works computed either single-site (or 
single block) von Neumann entropies or two-spin entan- 
glement measures such as concurrence Q and negativ- 
ity [1], and demonstrated that these measures or their 
derivatives display a peak (discontinuous jump) either 
near or at the second order (first order) quantum critical 
point (QCP). They can therefore serve as tools for iden- 
tifying quantum phase transitions in equilibrium quan- 
tum critical systems. These studies provide a bridge 
between quantum information theory and equilibrium 
quantum critical phenomena which can be useful for sev- 
eral aspects of quantum computations, cryptography and 
teleportation 9]. However, in all of these studies, the 
computed entanglement receives contribution from the 
ground state of the systems alone; excited states are not 
probed. 

Recently, enormous progress has been made in un- 
derstanding defect production due to non-equilibrium 
dynamics of a system passing through a critical point 
[iaiililil- In particular, it was shown that for a 
slow linear quench through a QCP, the defect density 
n scales with the quench time r with an universal expo- 
nent: n ~ ^-yd/(zv+i) ^ where v and z are the correla- 
tion length and dynamical critical exponents associated 
with the phase transition and d is the system dimension 
[Tllfl3] . Such results have been extended to cases where 
the quench takes the system through a quantum critical 
surface along a gapless line or a multicritical point 
[Tst . and for non-linear power-law quenches ^1^]. More 
recently, the two-spin entanglement properties of a quan- 
tum system during time evolution after a sudden quench 



through a critical point have been studied [iTj. However, 
the nature of two-spin entanglement generation due to 
a finite rate of quench through a quantum critical point 
has not been investigated so far, although other kinds of 
entanglement have been studied 18]. 

In this work, we compute concurrence and negativ- 
ity as measures of two-spin entanglement generated by a 
power-law quench, characterized by a rate r and an expo- 
nent a, which takes a spin-1/2 XY chain in a transverse 
field through a QCP. Our central results are as follows. 
First, we show that, in contrast to the studies of entangle- 
ment in anisotropic XY chains so far , such a quench 
generates entanglement only between spins separated by 
an even number of lattice spacings (even neighbors) ; the 
nearest neighbor sites are not entangled. Second, a crit- 
ical quench rate is required to generate two-spin en- 
tanglement (unlike other kinds of entanglement which do 
not seem to require a critical quench rate [Ts] ) . For faster 
quench rates (r < Tc), there is no entanglement between 
any pair of sites, and the entire entanglement is multipar- 
tite which is rather unusual. Third, by tuning the quench 
rate, one can control the amount of entanglement, and 
tune the ratio of entanglements between a spin and its 
consecutive even neighbors to take values between zero 
and 1 . Finally, for large quench time r and a given power 
a, the concurrence (negativity) scales as yja/r (a/r); 
this scaling is directly related to defect production by 
the quench through the QCP. To the best of our knowl- 
edge, the scaling behavior of the two-spin entanglement 
generated by a slow quench through a QCP and its selec- 
tive generation by tuning the quench rate have not been 
reported so far. Hence our study constitutes a significant 
extension of our current understanding of the nature of 
entanglement in many-body systems. This study may be 
relevant for the generation of entanglement as a resource 
for many algorithms of quantum computations, cryptog- 
raphy and teleportation Q . 

The organization of the rest of this paper is as follows. 
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In Sec. |TT1 we compute concurrence and negativity as 
measures of entanglement for the XY model. This is 
followed by a discussion of our main results in Sec. IIIII 
We end with some concluding remarks in Sec. IIVI 



II. MEASURES OF ENTANGLEMENT 

We begin with the Hamiltonian of the spin- 1/2 
anisotropic XY spin chain given by 

^ = I E[(i + + (1 - + h<], (1) 



where cr" for a = x,y,z are the Pauli matrices, J(l -I- 
[—)"/) /2 are interaction strengths between x{y) com- 
ponents of the nearest-neighbor spins (we set the lat- 
tice spacing d = 1), h is the magnetic field in units 
of J [we henceforth set h — J = 1, and measure 
all energies (times) in units of J (h/J)], and 7 is the 
anisotropy parameter which varies between (isotropic 
XY chain) and 1 (Ising limit). The phase diagram of 
the above model can easily be obtained by mapping the 
spins to fermions via a Jordan- Wigner transformation: 

a! = - 1/2, and cM) = ItP-c. '^j (-'^Y'^i i^t) 
[l9ll20| . The fermions can be shown to have an effec- 
tive two- level Hamiltonian for each pair of momenta ±fc 
in terms of the states |0) and \k,—k) — c[.cLj,|0) 
by iJfc = — [h + 2cosfc] (/ + r^) — 2Tj,7sinfc 



given 



where /, and Ty denote identity and Pauli matrices in 
the |0), \k,—k) space. The equilibrium phase diagram 
for the model is well-known; for \h\ > 1, there is a para- 
magnetic phase with (erf) = 0, while for \h\ < 1, one 
finds a ferromagnetic phase with (trf ) ^0. At h — ±1, 
there is a second order quantum phase transition with 
z = v = 1. At the transition dA, h — 1(— 1), the fermionic 
modes are gapless for k = 0(7r). The quench dynamics 
of the model, for a power-law time variation of the mag- 
netic field h{t) = /i|t/r|"sgn(t), where sgn is the signum 
function, has also been studied p^l20| . The quench starts 
with all spins down a,tt ^ —00 (i.e., the state |0) for all k) 
and ends ait—foo with a state in which the probabilities 
of |0) and |fc, — fc) are given by pk and 1 ~pk respectively. 
Here pk is the defect formation probability and is given 
by Pk = exp(-7rreff7^ sin^ k), where TeS = r/a [ll[2fl|. 

Armed with these results, we now compute the concur- 
rence and negativity as measures of the two-spin entan- 
glement of the spin chain generated by the quench. We 
note at the outset that the ground states of the initial 
and final Hamiltonians, at the beginning and end of such 
a quench process, are paramagnetic and do not posses 
any two-spin entanglement. Thus we expect that both 
for very fast (r — > 0) and very slow (r — > 00) quenches, 
where the system retains information only about the ini- 
tial and final ground states, the two-spin entanglement 
will vanish. Hence any finite entanglement obtained af- 
ter such a quench with a finite rate r must be generated 
by the non-adiabatic quench process and must therefore 



have contributions from excited states of the system. To 
compute the concurrence and negativity, we first note 
that the two-spin density matrix of the spin chain for 
any two sites i and j = i + n is given by |2l| 












^0 



"2 



a1 




(2) 



where the matrix elements a^., Oq and 6" 2 can be ex- 
pressed in terms of the two-spin correlation functions 

«± = (i(l±a,^)(l±<+J), 

a- = (i(l±<)(lT<+J), 



"1(2) 



-( + ) 



(3) 



The symmetry under af ^ — erf , af ^ —a^, 
ensures that all correlation functions such as 
and hence the remaining matrix elements are zero. The 
non-zero correlation functions for an arbitrary non-linear 
quench can be computed by generalizing the method de- 
veloped for a linear quench in Ref. [2^. We define 



an = — pkCOs{nk), 
Jo I" 



(4) 



and note that since pk is invariant under k —^ ir — k, 
a„ = when n is odd. In terms of a„, the diagonal 
correlation functions are given by 



1 



cr, a. 



z\2 



2aQ, 



(5) 



Thus, for any two spins separated by odd number of lat- 
tice spacings, {(^i crf_^_j^) — {crf)^ [20]. The off-diagonal 
correlators (cf c^+n) (where a, b can take the values 
— ) can also be computed in terms of a„. We shall present 
explicit expressions for these for n < 6 and provide a 
qualitative discussion for large n later. We find that 
(CTi^cTj^^) = &" = for all n since these involve cor- 
relations between two fermionic annihilation or creation 
operators and hence vanish. Further, {(^^o-f_^_^) = 63 = 
for all odd n since these are odd under the transforma- 
tion af ^ (-l)V^, ay ^ (-l)VK, < ^ which 



changes Jx.y — 
n < 6, we find 



— Jx,y and leaves pk invariant. For even 



(a4(a,f)-2a^)(af<+2), 
"6(«+2> + 4a2(a^ +al- a4«)) 
' ''"(TfcTf+2)-4(a2+a2)] + 16a^a4 



(6) 



Using these correlation functions, we can find all the non- 
zero matrix elements of for n < 6 and hence compute 
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FIG. 1: (Color online) Plots of C" as a function of r for n — 2 
(black solid line), n = 4 (red dashed line) and n — 6 (blue 
dotted line) and for 7 = = 1. The inset shows analogous 
plots for TV". 



FIG. 2: Variation of the position, r^, of the maxima of A/"^ 
as a function of a. The inset shows the values of the maxima, 
A/"™, as a function of a. 



the concurrence and negativity. The concurrence is given 
by C" = max{0, ^/X^ - ^/>^ - - where A^'s 

are the eigenvalues of p^^{ay p^* ^a^) in decreasing 
order Q. The \/Af are given by y^a"a" (appearing 
twice), and Oq ± |&2 1- Thus the spin chain has a non-zero 
concurrence if | &2 I > ^a" a" given by 

C" = max{0,2(|6«|-y^)} (7) 

To compute the negativity A/"", we need to take a partial 
transpose of with respect to the labels corresponding 
to the site j = i + n in Eq. ^ [8|. This interchanges 
^ 62 ; the eigenvalues of the resultant matrix p" are 
given byAg — (appearing twice) , and X^. = (l/2)[a" + 
a" ± ^ (a" — a" )^ + 4|62 P] of which only A" can become 
negative. This happens when |fe2 I > ^a"a" and yields 

AA" = maxIclA!!!} . (8) 

In the next section, we shall discuss the implications of 
these measures of entanglement in the context of XY 
model in a transverse field. 



III. RESULTS 

Eqs. ([7]) and ^ are the central results of this work 
which lead to several conclusions about the two-spin en- 
tanglement generated due to the quench for XY model 
in a transverse field. First, for odd n, (af tT|) = (cf )^ and 
(cr+cr") = 0. Thus all eigenvalues of /9"(cr2' ® aVp'^*ay (g) 
tj*') are equal to ao(l — ao) leading to C" = 0. All the 
eigenvalues of p are also positive; hence A/"" = 0. Thus 
the quench generates entanglement only between the even 



neighbor sites. Second, for large Toff and n <ti ^Tofi, 
a„ scales similarly to the defect density: a„ ~ ^JolJt. 
Using this, we find from Eqs. dSHS]) that 63 ^ \/a/T 
and a"a" ^ a/r which leads to C" ~ ^a/r and 
A/"" ~ a/r. Thus, for slow quenches, C" (A/"") scales with 
the same (twice the) universal exponent as the defect 
density [iJIilli^. This result relates two-spin entangle- 
ment generation for slow quenches to defect production 
by such a process. Third, both the concurrence and the 
negativity become non-zero for a finite critical quench 
rate (t")~^ (which is the solution of = a" a") 

above which there is no entanglement between a site and 
its n**^ neighbor. Solving this equation numerically, we 
find 7^^" = 1.96, 13.6 and 33.8, for n = 2, 4 and 6 when 
Q = 1. We shall shortly see that for r < r^, the entan- 
glement is entirely multipartite. 

We now plot C" and A/"" as a function of t for n = 
2,4,6 in Fig. [T] for 7 = a = 1. From Fig. [U we find 
that C" and A/"" becomes non-zero between r = r" and 
00. Further, the ratios C'/C^ oi N^/N'^ can be selectively 
tuned between zero and 1 by tuning r. The maximum val- 
ues of both C" and A/"" decrease rapidly with n. For large 
n ^ y^, using the properties of Toeplitz determinants 
used to compute the spin-correlators in these systems, it 
can be shown that {<7^cr~^^) ~ exp(— nZ-y/r) [201 • Thus 
we expect the entanglement to vanish exponentially for 
n ^ -y/r. In Fig. [21 we plot the positions and the mag- 
nitudes Af^ of the peaks in A/"^ as a function of a. We 
find that varies linearly with a while Af^ are indepen- 
dent of a. This behavior can be understood by noting 
that C" and A/"" depend on a through Teff. The peak 
of C" (A/") occurs when dC/dr = a''^ dC' / dr^s = 
[dN'^/dr = a-^dW/dTeS = 0). Thus the position of 
the maxima, which is a solution to this equation, is given 
by T^gfj = Qf/r^ and hence varies linearly with a. The 
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FIG. 3: (Color online) Plot of scaled negativity jV^{t)T as a 
function of tj ^ for a = 7 = 1 and several r = 14, 16, 18, 20, 
showing the time evolution of entanglement generated by the 
quench which started at ti — —250. 



FIG. 4: Plot of M as a function of r showing the evolution 
of the multipartite part of the entanglement with the quench 
time. 



maximum value of the entanglement V^m) ^ func- 
tion of T^cff a-lone; it does not change with a. 

Next, we study the time evolution of entanglement by 
computing A/"" at different times t during the quench for 
a = 7 = 1 and n — 2. The plots of the scaled negativity 
M'^{t)T as a function of the scaled time is shown in 

Fig. [3]for several t. From this plot, we find that for large 
enough t, the generation of the entanglement starts when 
the system is near the QCP and reaches an asymptotic 
value M"^ as the quench takes the system away from the 
quantum critical point. The collapse of all the asymp- 
totic values of the negativity to the same curve for all r 
clearly confirms the scaling behavior of the entanglement 
discussed earlier. Similar characteristics hold for A/"**'^. 

We now discuss the nature of the entanglement gen- 
erated by the quench. We note that the single-site den- 
sity matrix of the spin chain is pi = {I + a^{a^))/2. 
Using Eq. ([3]), we find that the single-site concurrence 
C(i) = V4detp, = 2y^ao(l - ao) [H is finite for all r ex- 
cept r = 0, 00 for which ao — 1,0. Since C^^-* is a measure 
of the entanglement of spin i with all the other spins in 
the chain, we conclude that the entanglement generated 
when r < is entirely multipartite. The multipartite 
par t of the entanglement for any r can be quantified as 



M 



(C 



(1)^2 



(9) 



n>2 



A plot of M as a function of r is shown in Fig. IH where 
we have summed C„ for n < 6. We find that M decreases 
rather slowly; hence a determination of the fate of M for 
large r appears to be a rather difficult task. We note 
that for large r, (C^i))^ ^ as l/^/r, while (C")^ ^ 
as 1/r for t ^ n^. Hence M may either go to as 



a power of r as t ^ 00, or may vanish at some large 
T beyond which all the entanglement becomes bipartite. 
Differentiating between these two possibilities is left as a 
subject of future study. 

An experimental verification of our theory would in- 
volve measuring the two-spin correlation functions of 
an anisotropic XY chain after performing a quench. 
There are several compounds such as K3Fe(CN)6 (J — 
-0.23K), (NH4)2MnF5 (J ~ -12K), and RbFeClg • 2H2O 
(J ~ — 35K) where the Ising limits of these chains are 
realized [2^. Similar experiments, involving measure- 
ment of two-spin correlation functions in equilibrium us- 
ing neutron scattering, have recently been carried out for 
square-lattice antiferromagnets, where short-range en- 
tanglement between spins has been demonstrated [24| . 



IV. CONCLUSIONS 

To conclude, we have shown that a controlled amount 
of entanglement can be generated in an anisotropic XY 
spin chain by performing a power-law quench; such a pro- 
cess generates two-spin entanglement only between sites 
which are even neighbors. The generated entanglement 
is entirely multipartite when the quench is faster than 
a critical rate; such states which have only multipartite 
entanglement are quite uncommon. The entanglement 
between even neighbors shows a scaling behavior for slow 
quenches, similar to the scaling of defects. We note here 
that a similar calculation can be performed for the two- 
dimensional Kitaev model 25| since the exact two-spin 
correlations functions of the model during a quench pro- 
cess with an arbitrary rate r have been computed in Ref. 

. Such a calculation shows that the Kitaev model has 
zero bipartite entanglement (C" = A/"" = for all n) for 
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all quench times t, and the entire entanglement is always 
multipartite. This leaves us with the question of the re- 
lation between bipartite entanglement generated during 
a quench and the integrability and number of conserved 
quantities of the underlying system; this would be an 



interesting subject for future studies. 
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